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We investigate the decay of B° — > K*°R and B l j —> <f>R with R being the X(4160), T(3940), 

Z(3930) resonances. Under the assumption that these states are dynamically generated from the 
vector-vector interaction, as has been concluded from several theoretical studies, we use a reaction 
mechanism of quark production at the elementary level, followed by hadronization of one final qq 
pair into two vectors and posterior final state interaction of this pair of vector mesons to produce 
the resonances. With this procedure we are able to predict five ratios for these decays, which are 
closely linked to the dynamical nature of these states, and also predict the order of magnitude of 
the branching ratios which we find of the order of 1CD 4 , well within the present measurable range. 

In order to further test the dynamical nature of these resonances we study the B° —» cpD* D* and 
B° —» (j)D*sDl decays close to the D*D* and D* S D* S thresholds and make predictions for the ratio of 
the mass distributions in these decays and the B'j —> f>R decay widths. The measurement of these 
decays rates can help unravel the nature of these resonances. 


I. INTRODUCTION 


The XYZ resonances with masses in the region 
around 4000 MeV have stirred the hadron community 
with a series of states that challenge the common wis¬ 
dom of mesons as made from qq. There has been in¬ 
tense experimental work done at the BABAR, BELLE, 
CLEO, BES and other collaborations, and many hopes 
are pinned in the role that the future FAIR facility with 
the PANDA collaboration and J-PARC will play in this 
field. There are early experimental reviews on the topic 
m and more recent ones EHS- From the theoretical 
point of view there has also been an intensive activity 
trying to understand these intriguing states. There are 
quark model pictures [10} [TT] and explicit tetraquark 
structures [T2]. Molecular interpretations are given in 
refs. HMD- Much progress also has been done using 
the heavy quark spin symmetry (HQSS) [211 [233 ■ P re " 
dictions using QCD sum rules have also brought some 
light into the issue pMH251 . Strong decays of these res¬ 
onances have been studied to learn about the nature of 
these states mm , while very often radiative decays 


are invoked as a tool to provide insight into this problem 
[2911331 , although there might be exceptions as discussed 
in ref. [55]. It has even been speculated that some 
states found near thresholds of two mesons could just 
be cusps, or threshold effects [35,. However, this spec¬ 
ulation was challenged in ref. [SB] which showed that 
the near threshold narrow structures cannot be simply 
explained by kinematical threshold cusps in the corre¬ 
sponding elastic channels but require the presence of 
.S-matrix poles. Along this latter point one should also 
mention a recent work calling the attention to possible 
effects of singularities on the opposite side of the uni¬ 
tary cut that enhance the cusp structure for states with 
mass above a threshold G3- Some theoretical reports 
on these issues can be found in refs. [551140] . 

On the other hand, and somewhat unexpected, recent 
experiments on the weak decays of B meson are proving 
to be a powerful source of information on hadron dy¬ 
namics and the nature of hadronic states [4TH46] . One 
of the recent surprises was to see from these experi¬ 
ments a pronounced peak for the /o(980) in decay 
into J/if and 7r + 7r - [41] while the signal for the /o(500) 
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was found very small or non-existent. Simultaneously, 
in the analogous decay of B° into J/ip and 7r + 7r“ [45] 
a signal was seen for the / o (500) production and only a 
very small fraction was observed for the /o(980) produc¬ 
tion. This is surprising since the /o(500) couples more 
strongly to 7r + 7r“ than to the /o(980). Some attempt to 
explain this behaviour was offered in ref. m in terms 
of tetraquark structures for the scalar mesons. Fur¬ 
thermore, the / 0 (500) and /o(980) resonances are nat¬ 
urally explained as states dynamically generated from 
the meson-meson interaction in chiral unitary theory 
g5H52|. From this perspective the features observed 
in these experiments and different ratios were well de¬ 
scribed in ref. [54]. The basic approach in ref. [54] was 
to identify the dominant mechanism of the decays at the 
quark level, implementing the hadronization of the final 
qq pair into two mesons and allowing them to interact 
to generate the /o(500) and / 0 ( 980) resonances. 

The picture of ref. [54] has been extended to de¬ 
scribe many other B and D decays, with equal success 
in the description of observed features of the reactions. 
In ref. [55] ratios for the production of vector mesons in 
the final states were evaluated and predictions for the 
B° s —>• J/ipn( 800) decay were made. In ref. [SB] the 
D° decays into K° and /o(500), /o(980), ao(980) were 
studied. A different sort of states, dynamically gen¬ 
erated from the vector-vector interaction were investi¬ 
gated in ref. m in the B° and B® decays into J/ip and 
/o(1370), /o(1710), / 2 (1270), /£(1525), K* 2 { 1430). In 
ref. [58] the B° decay into D° and p or / 0 (500), /o(980), 
ao(980) and decay into D° and K*° or k(800) were 
studied. Similarly, the decay into D s DK was stud¬ 
ied with the aim of learning about KD scattering and 
the D* 0 (2317) resonance [59] . Related work to this lat¬ 
ter one was done studying the semileptonic B s and B 
decays in ref. [BU]. 

Related problems of these weak decays have been ad¬ 
dressed in refs. [BMZI from a very different point of 
view, evaluating microscopically the weak matrix ele¬ 
ments and parametrizing parts of this interaction plus 
properties of the resonances produced, then carrying fits 
to data. The aims are also usually different, having in 
mind the construction of full amplitudes that can be 
used for issues like CP violation. On the contrary, the 
work of ref. [54] and the related ones make predictions 
for shapes and ratios that are tied to the final state in¬ 
teraction of mesons and, consequently, parameter free 
predictions can be made for these observables. 


So far, in the study of these B decays the production 
of XYZ states has not yet been addressed and the aim 
of the present paper is to study reactions where these 
states can be produced, evaluating ratios for different 
decay modes and estimating the absolute rates. This 
should stimulate experimental work that can shed light 
on the nature of some of these controversial states. 


II. FORMALISM 

Following refs. HUES, we plot in Fig. [T] the basic 
mechanism at the quark level for B®(B°) decay into 
a final cc and another qq pair. In ref. [51] the cc 



FIG. 1. Diagrams at the quark level for B® (a) and B° (b) 
decays into cc and a qq pair. 

went into the production of a J/ij) and the ss or sd 
were hadronized to produce two mesons which were al¬ 
lowed to interact to produce some resonant states. Here, 
we shall follow a different strategy and allow the cc to 
hadronize into two vector mesons, while the ss and sd 
will make the 4> and K*° mesons respectively. Let us 
observe that, apart for the b —> c transition, most fa¬ 
vored for the decay, we have selected an s in the fi¬ 
nal state which makes the c —> s transition Cabibbo 
allowed. This choice magnifies the decay rate, which 
should then be of the same order of magnitude as the 
B° s —»• J/if)f 0 (980), which also had the same diagram of 
the quark level prior to the hadronization of the ss to 
produce two mesons, in this case KK that couples later 
to the / o (980). 
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The next step consists in introducing a new qq state 
with the quantum numbers of the vacuum, uu + dd + 
ss + cc, and see which combinations of mesons appear 
when added to cc. This is depicted in Fig. [2] 

An easy way to see which vector mesons are produced 
in the hadronization of cc is to introduce the qq matrix 


M = 


l uu ud us uc \ 
du dd ds dc 
su sd ss sc 
cu cd cs cc 


(u d s c ) . (1) 


Note that this matrix corresponds to the SU(4) vector 
matrix 


V = 


/ + 

v^2 ' y/2 

p~ 
K*~ 

D*o 


p + 

K *+ 

D *o \ 

P° I U) 
\/2 ' \/2 

K*° 

D*~ 

K*° 


D*~ 

D *+ 

D*+ 

•W 


( 2 ) 


Now we see that [68] 


M-M = | “ ] ( u d 



u a s c 


u d s c ) (uu + dd + ss + cc) 


= M(uu + dd + ss + cc). (3) 

Hence, we can write 

cc(uu + dd + ss + cc) = (M ■ M) 44 = ( V ■ V) 44 (4) 

and 


(V ■ V) 44 = D*°D*° + D*+D*~ + D*+D*~ + J/if>J/i(>. 

, (5) 

Note that we have produced an I = 0 combination, 
as it should be coming from cc and the strong inter¬ 
action hadronization, since we have the isospin doublet 
(D* + , — D*°), (D*°,D*~). The J/'i/'J/'i/' component is 
energetically forbidden and hence we can write 

(V ■ V) 44 —> V2(D*D*) I=0 + D* s +D* s ~. (6) 


We can reach the result of Eq. ([ 5 ]) in a simpler way 
neglecting the cccc component in Eq. Q. Since cc is an 
SU(3) singlet and so is the uu + dd + ss combination, 
then the D*D* combination of Eq. (|6|) is also an SU(3) 
singlet. Using the full formalism has the advantage that 
it already tells one the phase convention of the states, 
which should be taken consistently with the chiral study 
of the interaction based on the matrix of Eq. ([ 2 ]) . 

Following the philosophy of ref. [53[ we shall now 
let these vector mesons undergo interaction and we can 
connect with the work of ref. ClTj| . where using an ex¬ 
tension of the local hidden gauge approach SMZ3] some 
XYZ states were dynamically generated. Concretely, in 
ref. [69] four resonances were found, that we summa¬ 
rize in Table |T] together with the channel to which the 
resonance couples most strongly, and the experimen¬ 
tal state to which they were associated. In ref. [69] . 


Energy [MeV] 

I°[J PC ] 

Strongest Experimental 

channel state 

3943 - i7.4 

0+[0+ +] 

D*D* 

F(3940) [761 

3945 - iO 

o-[i+ ”] 

D*D* 

? Y P 

3922 - i26 

0+[2+ +] 

D*D* 

Z(3930) [77] 

4169 - i66 

0+[2+ +] 

DID,* 

A(4160) ES] 


TABLE I. States found in ref. |69 , the channel to which 
they couple most strongly, and the experimental states to 
which they are associated (see also refs. [Tim])- Y P is a 
predicted resonance. 

another state with 1 = 1 was found, but this one can¬ 
not be produced with the hadronization of cc. Some of 
these resonances have also been claimed to be of D*D* 
or D* S D* molecular nature in refs. [T5MST] using for it 
the Weinberg compositeness condition [82U84] and in 
refs. jSSMET] using QCD sum rules. Also using HQSS 
the same conclusions are reached in ref. 188] and with 
phenomenological potentials in ref. ,•*!)]. 

The final state interaction of the D*D* and D*D* 
proceeds diagrammatically as depicted in Fig. [3] Start¬ 
ing from Eq. ([6]) the analytical expression for the for¬ 
mation of the resonance R is given by 

&R) = Vp('fi‘9D*D*,RGD-D*+9D*D*,RG D *D*), 

, (7) 

where Gmm 1 is the loop function of the two intermediate 
meson propagators and gMM'.R is the coupling of the 
resonance to the MM' meson pair. 
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FIG. 2. Hadronization of the cc pair into two vector mesons for decay (a) and B° 


decay (b). 


The formalism for B° —>• K*°R runs parallel since 
the hadronization procedure is identical, coming from 
the cc, only the final state of qq is the K*° rather than 
the 4>. Hence, the matrix element is identical to the one 
of B ° — > 4>R, only the kinematics to different masses 
will change. 

There is one last point to consider which is the an¬ 
gular momentum conservation. For Jp = 0 + ,2 + , we 
have the transition 0“ —I J p 1~. Parity is not con¬ 
served but the angular momentum is. By choosing the 
lowest orbital momentum L , we see that L = 0 for 
J p = 1 + and L = 1 for J p = 0 + ,2 + . However, the 
dynamics will be different for J p = 0 + ,l + ,2 + . This 
means that we can relate B® —> y(3940)<)> with B° —► 
T(3940)A'* 0 , B° s -t Z(3930)<^ with B° -> Z(3930)X*°, 
B° s ->• X(4160)<)> with B° -» X(4160)A* 0 and B° s -t 
Yp(j) with B° —> YpK*°, but in addition we can re¬ 
late B° —► Z(3930)(f> with B® —> X(4160)()>, and the 
same for B° -> Z(3930)A'*° with B° -> X(4160)A*°. 
Hence in this latter case we also have a 2 + state for 
both resonances and the only difference between them 
is the different coupling to D*D* and D*D*, where the 
Z( 3930) couples mostly to D*D*, while the X(4160) 
couples mostly to D*D*. 

The partial decay width of these transitions is given 

by 


r* = 


8tt m 2 B o 




' p2L+l 

qb(K * 0 )’ 


( 8 ) 


which allows us to obtain the following ratios, where 
the different unknown constants Vp, that summarize the 


production amplitude at tree level, cancel in the ratios: 


Ri = 

1 B°^nt>RJ =0 

ba 

to 

III 

1 

bso_j, if .opj=o 

^ B°^K*orJ=i 

ba 

w 

III 

r bo-h/,rj=2 

R4 = 

^B0^d>R^= 2 

^ B°^K-°R’= 2 ’ 

Fb°^>k*°rJ= 2 


and 

_ V B^4,R(= 2 
As = Try-> 

1 B°^k/>RJ= 2 

where R J=0 , R J=l , R( =2 and R^ =2 are the F(3940), 
Yp, Z(3930) and A'(4160), respectively. 

III. RESULTS 

The couplings gMM' ,r and the loop functions Gmm' 
in Eq. 0. are taken from ref. [69| . where the dimen¬ 
sional regularization was used to deal with the diver¬ 
gence of Gmm 1 , fixing the regularization scale fj, = 1000 
MeV and the subtraction constant a = —2.07. How¬ 
ever, we have taken advantage to make corrections to 
the work of ref. |(JS) due to the findings of ref. [50] con¬ 
cerning heavy quark spin symmetry. It was found there 
that a factor tod* /mje has to be implemented in the 
hidden gauge coupling g = m p /2f n in order to account 
for the D* —»• Dir decay. However, this factor should 
not be implemented in the Weinberg-Tomozawa terms 
(coming from exchange of vector mesons) because these 
terms automatically implement this factor in the ver¬ 
tices of vector exchange. However in the work of ref. [55] 
the coupling needed for D* —>■ Dir had been used also 
in the Weinberg-Tomozawa term. The bindings were 
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K*° 


(a) 


(b) 


FIG. 3. Diagrammatic representation of the formation of the resonances R(X, Y, Z ) through rescattering of MM' (D*D* or 
D^Ds) and coupling to the resonance. 


obtained by fitting subtraction constants in the G func¬ 
tion. We also do the same now, but with the present 
reduced interaction the G function becomes more nega¬ 
tive. We use now [i = 1000 MeV and a = —2.19. 

In Eqs. ([9])-( 10) we summarize the results that we 
obtain, 


i?i = 0.95, R 2 = 0.96, R 3 = 0.95, R 4 = 0.83, (9) 


and 


R 5 = 0.84. (10) 

As we can see, all the ratios are of the order of unity, 
however, the simplicity of these results should not be 
confused with a triviality. The ratios close to unity for 
the (j) or K*° production are linked to the fact that the 
resonances are dynamically generated from D*D* and 
D*D *, which are produced by the hadronization of the 
cc pair. The ratio for the J p = 2 + is even more subtle 
since it is linked to the particular couplings of these res¬ 
onances to D*D* and D*D* : which are a consequence of 
the dynamics that generates these states. Actually, the 
ratios R 4 , R- 2 , R 3 , R 4 are based only on phase space 
and result from the elementary mechanisms of Fig [TJ 
We would get the same ratios as far as the resonances 
are cc based. Hence, even of these ratios do not prove 
the molecular nature of the resonances, they already 
provide valuable information telling us that they are cc 
based. 

The ratio R$ provides more information since it in¬ 
volves two independent resonances and it is not just a 
phase space ratio. If we take into account only phase 
space then R$ ss 4 instead of the value 0.84 that we 
obtain. 


As for the absolute rates, we can establish an analogy 
to the B° s —► J/ipfo(980) decay, since the amplitudes 
of the quark level prior to hadronization are identical 
and both processes require hadronization of a qq pair 
into resonances. Hence, we estimate the branching ratio 
for the production of these states of the order of 10 -4 
US]. Another estimate can be done starting from B® —► 
J/ip(t>, which has a branching ratio of 1.07 x 10” 3 and 
reducing this rate by one order of magnitude which is 
the reduction factor that is found for hadronization in 
ref. [55] . In both cases we find branching fractions of the 
order of 10~ 4 , which are an order of magnitude bigger 
than many rates of the order of 10~ 5 already catalogued 
in the PDG [75]. 

Given the fact that the ratios R\, R 2 , R 3 , R 4 ob¬ 
tained are not determining the molecular nature of the 
resonances, but only on the fact that they are cc based, 
we propose a complementary test in the next section. 


IV. COMPLEMENTARY TEST OF THE 
MOLECULAR NATURE OF THE RESONANCES 


In this section we propose a test that is linked to the 
molecular nature of the resonances. We study the decay 
Bg —>■ (j)D*D* or B ° —> <f>D*D* close to the D*D* and 
D*D* thresholds. 

Let us now look to the process B® —► <f>D*D* depicted 
in Fig. [4] The production matrix for this process will 
be given by 

=Vp(V 2 + V2Gi f(l_n) + G 2 f(2->-l)), 

( 11 ) 
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FIG. 4. Feynmann diagrams for the D*D* production in B° decays. 


where 1 and 2 stands for the D*D* and D*D* channels 
respectively. The differential cross section for produc¬ 
tion will be given by [54] 


dr 

dM inv 




W> (12) 


where is the </> momentum in the B° s rest frame and 
Pd* the D* momentum in the D*D* rest frame. By 
comparing this equation with Eq. ([S]) for the coalescence 
production of the resonance in B® —> </> R 1 we find 


Rr= 


M\ 1 dT 

Ptj>PD* r R dM inv 
Mr p/(M inv ) t(B°^<pD*D*) 

471-2 P4, L+1 (Mr) 


(13) 


where we have divided the ratio of widths by the phase 
space factor p^pn* and multiplied by Mji. to get a con¬ 
stant value at threshold and a dimensionless magnitude. 
We would apply this method for the three resonances 
that couple strongly to D*D* (see Table [ij). In the case 
of the resonance R 2 with J = 2 that couples mostly to 
the D*D* channel (see Table [ij) we look instead for the 
production of D*D* S . 

For the D*D* production we have 


t(B®—x/>D* D*) — Vp( 1 + V^Gi t( 1_>2) + G 2 t(2->2)), 

(14) 

and we use Eq. ( fl3| ) but with D*D* instead of D* D* 
in the final state. We have evaluated Eq. (13) using 
the scattering matrices obtained in ref. [BjJj modified as 
discussed above, and then Eq. © and Eq. ( [l4| . The 
results are shown in Fig. [5j 


We can see that the ratios are different for each case 
and have some structure. We observe that there is a fall 
down of the differential cross sections as a function of 
energy, as it would correspond to the tail of a resonance 
below threshold. We should also note that in the case 
of D*D*, we have produced the 1 = 0 combination. If 
instead, one component like D* + D*~ is observed, the 
rate should be multiplied by 1/2. In the case of D*D* 
there is a single component and the rate predicted is 
fine. 

In order to estimate uncertainties, we modify the 
model of ref. |T5] by changing a bit the value of a in 
the G function expressions such that the values of the 
masses of the states change 5 MeV up or down. The 
results are shown in Fig. [6j 

We find differences of the order of 15% in the D*D*, 
D* S D* S mass distributions, which we accept as systematic 
uncertainties. 

Concerning the uncertainties of the ratios, since 
Ri, i? 2 , R 3 , Ra come from phase space, there are 
no changes for these ratios when using the second 
model. The ratio R$ is modified and we find now 
R 5 = 0.84 ± 0.02. 

V. CONCLUSIONS 

We have investigated the decays of B° —► K*°R and 
B° s -)• <j)R with R being the X(4160), E(3940), Z{ 3930) 
and a predicted J = 1 resonances. These decays have 
not been yet investigated. We estimate them to have 
branching ratios of the order of 10~ 4 . We used a model 
in which these states are dynamically generated from 
the vector-vector interaction, most notably the D*D* 
and D*D* states. Within this model we could predict 
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M lnv (L>/D s *)[MeV] 


FIG. 5. Results of R [f' of Eq. |l3| as a function of Mi nv (D*D*) for the first three resonances of the Table]!] (left) and 
as a function of Mi nv (D*D*) (light) for the fourth resonance of the Table [i] 



FIG. 6. Results of R^ 


of Eq. (131 as a function of M lnv (D*D“) for spin= 0. 


five ratios for the production of these states, where un¬ 
known dynamics of the production processes cancels in 
the ratios. The procedure used a mechanism in which 
we have quark production at the elementary level, fol¬ 
lowed by hadronization of one final qq pair into two 


vectors and posterior final state interaction of this pair 
of vector mesons. This mechanism has been tested suc¬ 
cessfully in many other decays and given us confidence 
on the fairness of the predictions made. Some ratios 
only tell us that the resonance is build from cc with or 







without extra hadronization. In order to offer some ex¬ 
tra test for the molecular composition of these states, 
we have evaluated the invariant mass distributions for 
the B° s —y cj)D*D* and B® —> (j)D* s D* s close to threshold 
and have made predictions for these magnitudes rela¬ 
tive to the width for the production of the resonances. 
The experimental investigation of these decay modes, 
and comparison with the predictions made, would shed 
light on the nature of these resonances and we can only 
encourage the implementation of such experiments. 
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